Collective Excitations of a Two-Component Bose Condensate at Finite Temperature 
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We compare the collective modes for Bose-condensed systems with two degenerate components 
with and without spontaneous intercomponent coherence at finite temperature using the time- 
dependent Hartree-Fock approximation. We show that the interaction between the condensate and 
non-condensate in these two cases results in qualitatively different collective excitation spectra. We 
show that at zero temperature the single-particle excitations of the incoherent Bose condensate can 
be probed by intercomponent excitations. 
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Introduction. For a Bose-condensed system with two 
internal components, one of the fundamental questions 
is whether one can distinguish a system with sponta- 
neous intercomponent coherence from one without it . 
At zero temperature, the ground state energies calcu- 
lated with the static Gross-Pitaevskii energy functional 
are the same for both cases 0]. The collective ex- 
citations calculated with the linearized time-dependent 
Gross-Pitaevskii (TDGP) equation are also identical for 
both cases. For example, the collective modes for a ho- 
mogeneous two-component system are 



± 



yJ(u>*-u>tf + 16gl 2 p lP2 (e%) 



(1) 



for both cases, as obtained from the linearized TDGP 
or the Bogoliubov theory HH0], where £°- = 4^, 

Luf = e-i \A + 2gapij is the Bogoliubov mode for a single 
component system, gn is the i-th intracomponent contact 
interaction, gi% the intercomponent contact interaction, 
and pi is the density of i-th component with i — 1,2. The 
excitations in Eq. Q may be interpreted differently for 
the incoherent and coherent cases are different. Without 
intercomponent coherence, there are two macroscopically 
occupied states, a situation that has come to be known 
as a fragmented condensate jjj. In this case ui+ may be 
interpreted as a density wave in which the two conden- 
sates move in phase, while lo- is a collective excitation 
associated with out of phase motion between the two 
condensates. With intercomponent coherence, there is 
only a single condensate. In this case w + is interpreted 
as the superfluid mode, and oj_ is a Goldstone mode 
corresponding to a spontaneously broken U(l) symme- 
try in the relative phase between the two components. 
As pointed out by Leggett [lj, in the absence of parti- 
cle exchange between the two components, no physical 
quantities can depend on the relative phase between the 
them. However, at finite temperatures, we will show that 
the coupling between the condensate and non-condensate 
particles provides a way of probing the nature of the con- 
densates. 

In this work, we compare the collective excitations for 
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FIG. 1: (Color online) The appearing of the new out-of-phase 
mode in the incoherent two-component systems at finite tem- 
perature. 



a two-component Bose-condensed system with and with- 
out intercomponent coherence at finite temperature us- 
ing the time-dependent Hartree-Fock (TDHF) approxi- 
mation, for the special case <?n = 322 > 312- Our princi- 
pal findings are as follows. Without intercomponent co- 
herence: (i) The motions of the two components and un- 
condensed particles give rise to gapless collective excita- 
tions. We find one in-phase (superfluid) mode. However, 
there can be two out-of-phase modes within some range 
of intercomponent interaction strengths at appropriate 
temperatures. This is illustrated in Fig. ^ (h) The per- 
turbative response to a field that allows intercomponent 
conversion has a resonance at frequencies given by the 
single-particle excitation energies. With intercomponent 
coherence, we find a gapless superfluid mode, and two 
"pseudo-spin waves" , one gapless and the other gapped, 
as illustrated in Fig. [21 The gapped pseudo-spin wave 
has vanishing weight as temperature goes to zero. The 
differences in in these spectra allow one to distinquish a 
single condensate from a fragmented one. 

Model. We consider the collective excitations in a two- 
component homogenous Bose-condensed atomic gas at 
finite temperatures ksT = 1/(3. To simplify the physical 
picture, we assume the two internal states are degener- 
ate and take the two-body interactions as contact forms 



with strength parameters g\\ = 322 = 9s = 



and 
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FIG. 2: (Color online) The two pseudo-spin modes in the 
system with intercomponent coherence. 



.912 = 9x = 4 ^ a x with a s _ x the corresponding s-wave 
scattering length, and take g s > g x , so that the homoge- 
nous state is stable against phase separation. The system 
may be characterized by the gas parameter pa? s , where p 
is the total particle density, and 7 = g x /g s . Our unit 
of energy is Eq — ^Jna 1 ' wnere m the mass of the par- 
ticles. Using a pseudo-spin language, the single-particle 
wavefunctions may be written as spinors with compo- 
nents V'jf.x(r) = y ^e lk ' r rj\, with V the volume of the 
system and rj X the spin state functions. In this work, the 
spinor index is denoted by Greek letters. For incoher- 
ent condensates, we write A = 1,2 with 771 = 1 
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the pseudospin states of the two condensates. 
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be in the 77+ state. In both cases the macroscopically 
occupied states have zero momentum. 



For coherent case, we write A = ±, 77+ 
r/_ = -7- [771 — 772] , and the condensation is assumed to 



Method. We briefly outline our method here; de- 
tails will be presented elsewhere. It is well known that 
to treat the dynamics of a condensate and thermally 
excited (non-condensate) particles in a fully consistent 
manner is theoretically challenging ||. Recently we 
have shown this can be done within a TDHF approx- 
imation |9l by using a constrained grand canonical en- 
semble [lOj in which the condensate particle number is 
fixed. We first determine the static Hartree-Fock (HF) 
ground state properties, including the chemical poten- 
tial /1, the single-particle energies e^ Q , the occupation 
numbers n ? , and the critical temperature (3 C - This in- 
volves defining a self energy S which is a functional of 
the densities Pk lXl ,k 2 x 2 =< a l lXl a k 2 \ 2 >> with a kx the 
boson destruction operator, and the condensate wave- 
function^) t/>oa- We then introduce a small perturbation 
SU , and assume £ keeps its functional form with respect 
to the p's and ipo\, which are now time-dependent. This 
results in self-consistent equations which can be solved 
in the linear response regime. This allows us to ob- 



tain equations for the response functions X\a a/3 
"^(ot , (*Wx(i)4 (OW ,~(oj},' 



with a vertex 
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(3) 



Here s and s' denote the components of each spinor, and and one obtains 



^AmiA iA , = £ f i? ?A(s)?7 Ml (s')?7A 1 (s')?7 A1 (s). The Kronecker 
delta terms 5 Xc are 1 when A is a condensate state, and 
otherwise. Poles of the response functions correspond 
to collective excitations of the system. 

Results. We first consider incoherent conden- 
sates. In this case, the density-density response func- 
tion matrix has only six non-zero matrix elements: 

xii,ii(<T;w) = X22,22(<7*;^), xn,22(?;w) = X22,n(<f;w), 

and xi2,2i(q;w) = X2i,i2(-<f; -o>)> where XXn.a/siftw) = 

?Ek 1 ,fc 2 ZAfto/3(9»- Thc equations for the first four 
functions are decoupled from those of the latter two. 

The functions xn,ii(<f; and X22,n(s; w) describe the 
density response of the two components. From Eqs. @ 



1 - 2g s P -g x P \ f X u,u 
-g x P 1 - 2g s P I I X22,ii 



(4) 



where P(q;oS) = P c {q;w) + P n (q;uj) with P c (q;uj) = 
2 po4 nnr i p ( ff..,\ _ iv, 

V q l asru g fc + g fc 

with Tig = = ntn- The poles of these x's occur when 
the determinant of the matrix in Eq. @ vanishes, 

1 - (2g s ± fo ) [P c (g; w) + P n (q; u)] = 0. (5) 

Equation (JjjJ shows the effect of intercomponent inter- 
action and the interaction between the condensate and 
non-condensate on the collective excitations. The plus 
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and minus signs in Eq. JSJ) can be interpreted as result- 
ing from the in-phase and out-of-phase motion of the 
two components, respectively. We find only one in-phase 
(gapless) mode as a solution, which reduces to oj + in Eq. 
i[TJl in the zero temperature limit. Thus the interaction 
between the condensate and non-condensate has little ef- 
fect on the in-phase motion. However, this interaction 
has a striking effect on the out-of-phase motion. Within 
a range g Xj i < g x < g x ,2, the values g x ,i and g x ^ depend- 
ing on the gas parameter pa 3 and temperature T, the 
system can support a new out-of-phase mode for small q. 
This is illustrated Fig. □ for pa 3 = 1CT 4 and /3//3 c = 3, 
where it can be seen that |3 [xnn(9;w) — X22ii(<f; ^)] 
carries an extra resonance. In order to see why this oc- 
curs, we show the polarization functions P c and P n in 
Fig- El For £^ < %9sPo, P c {q\u) is a monotonic function 
of uj and approaches l/g s with zero slope as u> — > 0, while 
P n (q;w) is non-monotonic at small ui. The combination 
g s (2 — ~/)P(q]Lu) is then non-monotonic at small to, so 
for certain ranges of 7, Eq. 0is satisfied twice. This is 
shown in the inset in Fig. El 

The function X2i,i2(<f;w) supports a collective mode 
[ill Hzl Il3| , corresponding to an excitation in which the 
"flavor" of a particle changes. From Eq. J2Jl one finds 



X2i,i2(<f;w) 



where 
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2/30 [Ef + gsPO 



(6) 
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9sPo 



(7) 



At zero temperature the collective mode requires exciting 
a particle from one of the condensates to an excited state 
of the other component, so X2i,i2 has a pole at the single 
particle excitation spectrum, u> — e$ — p — e-i + g s po- 
Thus, this provides a way to probe the single-particle 
spectrum. We note that, since X2142 is the response to 
a "flavor-changing" perturbation, this collective mode is 
not found in the TDGP approximation when species ex- 
change is absent from the unperturbed Hamiltonian. It 
is apparent that our TDHF approach is able to capture 
a broader variety of excitations for this system. 

We now turn to the coherent case. The system can con- 
dense into any spinor given by rj(6) = ^= (771 Te^r^), 
resulting in a single condensate. We choose the system 



condensing into 77+ = rj(0 = 0), which spontaneously 
breaks the U(l) symmetry. Therefore, we would expect 
at least two gapless modes: a superfluid density mode 
and a Goldstone mode (a pseudo-spin wave) for the spon- 
taneous U(l) symmetry breaking. 

The density wave may be found from the poles of 

x++(q;oj) = Ek 1 k 2 ; af 3^+ 1 +%p(^ UJ ) and x— = 
1 
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FIG. 3: (Color online) The condensate and non-condensate 
polarization functions at which the new out-of-phase mode 
appears. The inset shows g s (2 — y)P(q;u>). The arrows show 
the locations of the two out-of-phase modes. 



^2k!k 2 -ai3 X k -- 2 a pi which follow a matrix equation ob- 
tained from Eq. J5J, 



1 - 9 P P++ -9sP+-i 
-g s P— 1 - g P P- 



P + 



, (8) 



where g p = g s + g x , P ++ {q;uj) = P++{q;uj) + 
P+ + (q;u>), P — (q; ui) = P?M") with P++(<t» = 

where the prime excludes the momentum of the conden- 
sate modes. As in the single component case, the non- 
condensate is too dilute to sustain a propagating second 
sound mode within this approximation. The resulting re- 
sponse and superfluid mode are thus qualitatively similar 
to the single component case. 

By contrast, the non-condensate component has a dra- 
matic effect on the pseudo-spin waves. In order to see 

this, we define X± T {q\u) = J2k 2a0 X±^ af3 (^ <*>) and 
split these into condensate and non-condensate parts: 

x+_(<f;w) = x+-°(<?;^), x-+(q;u) = x-+ q {q\u), and 

X±zp(<7;w) = Yl!k X±^(Qj w ) i where the prime on the sum 
excludes the momentum of the condensate mode. In 
terms of these quantities, Eq. (|3J| yields 
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1 - 9sPn g m P+-(q;ui) 
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^ + ~(<?» 


Vx+-(< 




\p-+(q;Lo)J 



(9) 



where g m = g s — g x . Eq. yields two spin waves, as 
shown in Fig. |3 one of which is gapless and originates 
from the U(l) symmetry breaking of the relative phase 
of the two internal states. At T = 0, this gapless mode is 
reduce to given by Eq. (JJJ. The other spin wave has 
a gap as q — > given by 

A = ^ [p§ + (p+ - p-) 2 ] + 9x9 p po(p+ - p-)- (10) 

Here p+ and /?_ are the non-condensate density in spinor 
r/ + and 7j_, respectively. The gapped spin wave, unlike 
the gapped flavor-changing mode obtained from Eq. © 
in the incoherent case, disappears as T — ► (i.e., its 
weight in the response functions vanishes), which can be 
see directly from Eq. © since both P"_ and P™ + are 
zero in this limit. 

In the above calculations we have discussed the special 
case of interaction strengths gn = #22 ■ For ju / 322 but 
.9i2 < 3ii322 so that the homogeneous state is still stable 
against phase separation, we find for the incoherent case 
results which are qualitatively the same. The coherent 
system, by contrast, is more complicated. In particular 
one finds that the spinors of the non-condensate modes 
are neither parallel nor anti-parallel to the condensate 
mode spinor. This leads to the response functions in 
Eq. all being coupled together, so a classification of 
the modes as flavor-changing and flavor-preserving is no 
longer possible. Nevertheless, one still finds two gapless 
linear modes and one gapped mode as we found in the 
3n = 322 case. Details of this more complicated situation 
will be presented elsewhere. 

In conclusion, we have compared the collective exci- 
tations of two-component Bose-condensed systems with 
and without coherence at finite temperature. The cou- 
pling between the condensate and non-condensate de- 
pends on whether there is intercomponent coherence, and 
the two identical dispersions given by Eq. JTJ at zero 



temperature evolve to quite different structures. We 
also have shown that, at zero temperature the single- 
particle excitations in the incoherent case can be probed 
by flavor-changing excitations, which can not be done in 
the coherent case. 
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